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Abstract 



We demonstrate the possibility of forming a single, large domain of dis- 
oriented chiral condensate (DCC) in a heavy-ion collision. In our scenario, 
rapid initial heating of the parton system provides a driving force for the chi- 
ral field, moving it away from the true vacuum and forcing it to go to the 
opposite point on the vacuum manifold. This converts the entire hot region 
into a single DCC domain. Subsequent rolling down of the chiral field to its 
true vacuum will then lead to emission of a large number of (approximately) 
coherent pions. The requirement of suppression of thermal fluctuations to 
maintain the (approximate) coherence of such a large DCC domain, favors 
three dimensional expansion of the plasma over the longitudinal expansion 
even at very early stages of evolution. This also constrains the maximum 
temperature of the system to lie within a window. We roughly estimate this 
window to be about 200 - 400 MeV. These results lead us to predict that 
extremely high energy collisions of very small nuclei (possibly hadrons) are 
better suited for observing signatures of a large DCC. Another possibility is 
to focus on peripheral collisions of heavy nuclei. 
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I. INTRODUCTION 



Recently there has been a lot of interest in exploring the very interesting possibility that 
extended regions, where the chiral field is misaligned from the true vacuum, may form in 
heavy- ion collisions or in large multiplicity hadronic collisions [0-^ • Such a region is termed 
as a disoriented chiral condensate (DCC). If a large DCC domain forms, then it can lead 
to spectacular signatures such as coherent emission of pions which can be detected p as 
anomalous fiuctuations in the ratio R of neutral pions to all pions. A strong motivation for 
looking for DCC comes from Centauro events in cosmic ray experiments 0. 

However, so far it has not been possible to realize this beautiful idea in the sense of 
finding experimental situations where large enough DCC domains can form which can lead 
to clear experimental signals. Essentially all theoretical estimates predict rather small DCC 
sizes, of the order of couple of fm at best. Not only is the number of coherent pions very small 
for small DCC, for a large QGP region undergoing phase transition, one expects many such 
DCC domains to form, making it harder to experimentally detect the signal from each DCC 
domain . Though, many clever ideas have been proposed to overcome these problems , 
one is still far away from realizing the ideal situation of having a large single DCC domain 
with clean experimental signature. 

In this paper, we demonstrate the possibility of forming such a large, single, DCC domain. 
In our model, the entire region of the plasma gets converted into a single DCC domain, which 
can have diameter as large as 14 fm for gold nucleus. We work within the framework of 
linear sigma model, and study the dynamics of the chiral field using Langevin equations of 
motion to take into account the effects of thermal fluctuations. The main idea of our model 
is that during the rapid heating of the initial parton system, the chiral fleld is driven towards 
the value zero due to a rapidly changing effective potential. If heating is fast enough, and 
the maximum temperature reached is large enough, then the fleld picks up enough kinetic 
energy to overshoot the zero of the fleld and goes to the opposite point on the vacuum 
manifold. That is, starting from the true vacuum with chiral angle being zero, the fleld ends 
up at chiral angle equal to tt. Most importantly, this happens for the entire region being 
heated up. The end result being that the entire region of the QGP ends up with chiral fleld 
being maximally disoriented, with the chiral angle being (close to) vr. Eventually the fleld 
will roll down from there, emitting coherent pions in the process. 

Due to thermal fluctuations, there will be some variation in the chiral angle over the entire 
DCC. The requirement that this variation remains relatively small, compared to the average 
disorientation of the chiral fleld from the true vacuum, leads to constraints on the maximum 
temperature, as well as on the rate of expansion of the plasma. If the maximum temperature 
To of the plasma is too large, or the rate of expansion too slow (e.g. in the longitudinal 
expansion model), then the system spends signiflcant amount of time at high temperatures 
so that thermal fluctuations become dominant. We flnd that thermal fluctuations remain 
small only when spherical expansion is assumed even at very early stages of the plasma 
evolution. For that case we flnd allowed values of Tq to lie roughly in the range 200 - 400 
MeV, the lower bound arising due to the fact that for small Tq the chiral fleld does not 
overshoot at all. Such values of Tq are small compared to the maximum QGP temperatures 
expected in central events in the heavy-ion collisions involving very heavy nuclei at LHC and 
RHIC. This suggests that in such experiments, or for that matter, even in present heavy- 
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ion experiments, observing a large DCC is more likely if one explores a range of centrality 
in selecting events. This is because with varying centrality the energy density, and hence 
plasma temperature, will also be expected to vary. Then the possibility will be larger for the 
plasma temperature to lie in the allowed range for some events. Requirement of early three 
dimensional expansion suggests that possibility of a large single DCC domain may be larger 
for very energetic collisions of smaller nuclei, or even hadron-hadron collisions. In this sense, 
Tevatron may indeed be a good place to look for DCC j^. This is, in some sense, consistent 
with what one expects from Centauro events where very heavy nuclei are not involved in 
the collision. Precise constraints on the masses of colliding nuclei and their energies, in the 
context of our model, are difficult to estimate at this stage. We hope to give somewhat more 
definite numbers in a future work. 

The paper is organized in the following manner. Sec. II discusses the physical motivations 
and various assumptions underlying the basic picture of our model. Sec. Ill discusses the 
numerical calculations and results. Our results naturally lead to consideration of the effects 
of the topology of field configurations on observables like pion distribution. This is discussed 
in Sec. IV. Conclusions are presented in Sec. V where we summarize our results and discuss 
various predictions. 



II. PHYSICAL PICTURE OF THE MODEL 

In this section we will discuss the basic picture underlying our model. Our numerical 
technique will be similar to that used in ref. and we will describe it in Sec. III. We 
will work within the framework of linear sigma model, with the finite temperature effective 
potential at one loop order given by, 

where represents components of the chiral field with $j = tTj for i = 1,2,3 and $4 = 
a. Here T is the temperature. We take = 140 MeV, = 94.5 MeV, and X 55 
(corresponding to rrio- = 1 GeV). In the chiral limit there is a second order phase transition 
with the critical temperature Tc = V^fw (— 134 MeV for our choice of parameters). With 
non-zero quark masses, chiral symmetry is explicitly broken, and is only approximately 
restored at a temperature ~ 115 MeV in the sense that the saddle point at the chiral angle 
6 = IT disappears above this temperature [1^. In the chiral limit, spontaneous breaking of 



chiral symmetry for temperatures below implies that one particular point on the vacuum 
manifold S*^ (characterized by J2i = fn) ^^^^ be chosen as the vacuum state in a given 
region of space, with all points on being equally likely. One may expect this to essentially 
hold true even in the presence of small pion mass. This will lead to a sort of domain structure 
in the physical space where each domain will have the chiral field aligned in a given direction, 
but the directions in different domains vary randomly. 

This essentially summarizes the conventional picture of the formation of DCC domains. 
Formation of this type of domain structure in a phase transition has been extensively dis- 
cussed in the context of topological defects in condensed matter physics and also in particle 
physics models of the early Universe. For an equilibrium, second order phase transition one 
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expects maximum size of domains to be of the order of m~^. For such small domains, 
dramatic signals like fluctuation in the ratio R of neutral pions to all pions will not be 
observable. However, it has been suggested that in a non-equilibrium transition coherent 
pion emission may be observable 0, see also ref. in this context. The amplification of 
long-wavelength pion modes due to parametric resonance, leading to larger DCC domains. 



has been discussed in ref. |12|. The effects of fluctuations on the growth of DCC domains. 



in the conventional picture, have been extensively studied in refs. IT^JT^. More recently. 



similar investigations, for the case of spherical expansion of the plasma, have been carried 
out in ref. ||T^ . Especially in the annealing scenario it is suggested that reasonably large 
DCC domains of sizes as large as 5 fm may form. The issue of initial conditions for these 
scenarios, appropriate at high temperatures, has been discussed in refs. Jl^, see also ref. 
T8| in this context. 



The most important aspect of all these models is that one starts with an intermediate 
chirally symmetric state, at temperature above T^. As the plasma cools due to expansion, 
and chiral symmetry breaks spontaneously, one considers small correlation domains and 
studies their growth. This inevitably leads to multi domain structure of the entire QGP 
region. 

We propose to choose a different starting point. Just after the collision of the nuclei or, 
say, just before it, the chiral field in the region is definitely in the true vacuum with the 
chiral angle equal to zero, and is described by the zero temperature effective potential. As 
the partons produced during the collision thermalize, eventually leading to a hot region with 
a temperature Tq (which will be in the chirally symmetric phase for large Tq), the effective 
potential changes to the one appropriate for that temperature. What happens during this 
time? It seems very hard to address this question, though, one can adopt different physical 
pictures modeling this intermediate regime. Presumably the most natural approach is to 
think of the initial system of partons as being completely out of equilibrium which then 
thermalizes due to interactions |T^J2^. As thermalization proceeds, the system approaches 
a state in (quasi) equilibrium with a well defined temperature, which is the maximum tem- 
perature of the QGP. Thus, in this picture, at any instant during this intermediate stage, 
the system does not have a well defined temperature. However, departure from equilibrium 
decreases, in the sense of various distributions gradually approaching the one appropriate 
for a system in thermal equilibrium at temperature Tq, as partons interact with each other 
and thermalize. 

Unfortunately, it is hard to do any calculations within this picture, as one has to deal with 
a non-equilibrium system, which is approaching an equilibrium state. Thus, for intermediate 
stages, use of finite temperature effective potential etc. becomes meaningless. Since the 
physics of our model is best described in terms of an evolving effective potential, we do 
not follow the above approach. Rather, we assume a simplified picture wherein we think 
of the initial system of partons at zero temperature which is getting heated up, finally 
reaching the temperature Tq of the QGP. In this picture, at any instant, one has a well 
defined temperature which first increases to a maximum value of Tq, and thereafter starts 
decreasing due to continued plasma expansion. The second part, that is cooling due to 
plasma expansion, is common to this picture as well as the conventional picture. However, 
initial heating assumed in the latter picture may appear confusing as the parton system is 
still expanding (in fact, rather fast in the beginning). Initial heating is actually supposed 
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to represent the thermalization of the parton system. So, even though plasma is expanding 
at any time, rate of thermahzation as well as secondary parton production dominates in the 
beginning, and leads to a consistent picture that the temperature is rising in the beginning. 
Once thermalization is complete, then only effect of expansion remains, leading to decrease 
in the temperature. Another way to think of the initial heating is to realize that we are only 
interested in the dynamics of the chiral field. It is clear that the chiral Lagrangian as given 
in Eqn.(l) does not account for all the degrees of freedom, especially at high temperatures. 
Thus, one can think of the interaction of these other degrees of freedom with the chiral field 
as giving rise to the heating of the system of the chiral field. 

To summarize the picture so far, we think of the initial thermalization of the parton 
system as, the system getting heated up to a maximum temperature Tq which subsequently 
cools due to continued plasma expansion. This translates in terms of the chiral model with 
a finite temperature effective potential with temperature being prescribed as a function of 
(proper) time. One can object on using the chiral model at very early times [0], and use of 
finite temperature effective potential for very short time scales in the pre-equilibrium phase 
of the parton system. However, we believe that the basic physics of our model has generic 
features which do not depend on such approximations. The spirit of our calculations is that 
by making these assumptions, we are able to demonstrate the essential physics of a novel 
mechanism for the conversion of the entire QGP region into a single DCC 

In order to give a clear picture of the physics underlying our model of DCC formation, 
we first consider a simple case with zero pion mass in the absence of any damping (or 
thermal fluctuations). Assume that the system initially is at zero temperature with the 
chiral field pointing in some direction on the vacuum manifold denoted by point P in Fig.l, 
all such directions having same energy for zero pion mass. Thick solid curve in Fig.l denotes 
the effective potential (in Eqn.(l)) at this stage. Now suppose that the temperature is 
instantaneously increased to a value Ti and kept there, where Ti is less than the critical 
temperature Tc. [For our choice of parameters, Tc ~ 134 MeV for zero pion mass.] The 
dashed curve in Fig.l denotes the effective potential at that temperature. As the temperature 
has been increased suddenly, the chiral field will still be having the value shown by point P 
in Fig.l. However, now the effective potential being given by the dashed curve, the chiral 
field at point P becomes unstable and starts rolling down towards the minimum of the new 
effective potential. It is clear that if the value of the potential energy at point P, with respect 
to the new effective potential, is larger than the height of the central bump then the chiral 
field will overshoot the bump and roll down to the opposite side on the vacuum manifold. 
From Eqn.(l) it is easy to see that this happens when Ti > Tc/V2. This then disorients 
the entire region of the plasma where the field is undergoing this dynamics. Of course, in 
this simplified model the chiral field will start oscillating about zero. In realistic case, due 
to cooling of the plasma, as well as due to damping from the expansion of the plasma (and 
from thermal fluctuations, particle production etc.), the field will get trapped on one side of 
the bump, with both sides being equally likely if several oscillations happen before the field 
settles down. We will show below that this basic picture works in rather generic situations, 
even with non-zero pion mass, and leads to conversion of the entire region (the region which 
has the same temperature evolution as a function of proper time) into a single DCC, apart 
from the effects of thermal fluctuations which we will discuss below. 
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III. NUMERICAL CALCULATIONS AND RESULTS 



We now discuss the realistic case with non-zero pion mass and with a temperature evo- 
lution consisting of initial heating stage with subsequent cooling. It is not clear what sort of 
time dependence is appropriate for the temperature in the thermalization (heating) stage. 
What we will do is carry out calculations for different types of time dependencies for the tem- 
perature during the heating stage, and present results for all such cases. For the temperature 



evolution during the cooling stage, we take the following time dependence [21 



T(r)=To(^)i. (2) 
r 

Here, Tq is the proper time at which the maximum QGP temperature Tq is achieved. We 
consider the longitudinal expansion model of Bjorken |2T[] for the plasma for which t] = 1, as 



well as the spherical expansion case where r] = 3. What we find is that spherical expansion 
is preferred for our model due to requirement of suppression of thermal fluctuations. Fig. 2 
shows plots of temperature as a function of the proper time for the spherical expansion 
case. We have chosen widely different possibilities for the temperature evolution during the 
heating stage, starting from the linear heating, shown by the dotted curve, to instantaneous 
heating with two different values of tq shown by the solid, and dashed curves. For the linear 
heating case, we keep the temperature equal to zero up to a value of r = rj„ which sets the 
scale for the expansion of the volume, as discussed below, = 0.2 fm and tq = 0.4 fm for 
the dotted curve in Fig.2.) 

To account for the effects of thermal fluctuations on the evolution of the chiral field, we 
study the dynamics using the following Langevin equations of motion 

+ + V')^^ = - ^ + Ur). (3) 

Here V{t) is the effective potential in Eqn (1). rj' is the friction coefficient due to coupling 
to the heat bath and ^i(T) represents Gaussian white noise term with < ^^(t) >= and 
< ^^(^2) >= ^v'^ij^i^i ~ ^2)- Here, V is the volume of the region over which the 

chiral field is being evolved by above equations. As mentioned in ref. [|],^, these Langevin 
equations can be interpreted as representing the dynamics of the soft modes (say zero mode 
over volume V) interacting with thermal fluctuations (heat bath) obtained by integrating out 
the hard modes which are assumed to be in thermal equilibrium, see ref. p2|j8[1 for details. 
In the absence of any thermal fluctuations, the chiral field, everywhere in the entire plasma 
region, Q, will end up at = tt after overshooting. However, in the presence of thermal 
fluctuations, the chiral field in different portions in Q, starting everywhere with 6 = 0, will 
overshoot to different points near 6 = it, thereby forming a patch at ^ = tt. We would like 
to get an estimate of the variation of the chiral field in Q, after overshooting, arising from 
thermal fluctuations. For this, we think of the whole plasma region Q in terms of smaller 
volumes V, and study evolution of average field in volume V with different realizations of 
the noise term in Eqn. (3). Each such realization then represents the evolution of the chiral 
field in one of the volumes V in the region Q. It is the issue of how the system becomes 
uncorrelated due to thermal fluctuations during a non-equilibrium evolution. There is a 
similarity between this and the situation studied in ref. where growth of domains was 
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studied for a system undergoing a quencli. In our case, tlie situation is reversed, as liere one 
wants to know liow correlations break up as the system is heated up. 

For the friction coefficient r^', an estimate from ref. p2|,p| for our case (with A ~ 55) will 



suggest a value of about 40 fm~^ (corresponding to the largest value of the temperature 
which is 350 MeV in our case , r/' will be smaller for lower temperatures). This value is far 
too large and the field dynamics is extremely damped in this case, and completely dominated 
by thermal fluctuations. For the parameters used in ref. the estimated value was about 
2 fm~^. As mentioned in ref. P], such large value of r]' gave large thermal fluctuations, 
hence small values (by about a factor 1/4) were also used in that work which gave smaller 
thermal fluctuations. In our case when we use a value of r^' ~ 4 fm~^, which is about 1/10 of 
the estimated value, we get decent results, though still thermal fluctuations are somewhat 
large. As we will argue below, in such case formation of large DCC is somewhat less likely, 
though our scenario of large DCC formation still works in this case. We find that thermal 
fiuctuations are under control when Vj is taken to be about 1 fm~^ or less. In this case our 
scenario works very well, and leads to formation of reasonably coherent, large DCC domain. 
Also, the probability of formation of such DCC does not suffer from extra suppression due 
to large thermal fluctuations. This suppression is really small for r/' = 0.6 fm~^ and, in this 
paper, we give detailed plots etc. for this value. We will quote results for the case when r]' 
is increased to larger values (e.g. 1, 4 and 10 fm~^). 

We now argue that the estimates of ref. p2|,|8| should not be rigidly applied for the 
strongly coupled case of linear sigma model and it is reasonable to consider much smaller 
values of r( This is because the calculations in ref. apply for a weakly coupled theory 
where perturbation calculations can be trusted. This is far from true in the present case 
where the coupling constant is as large as 55. (Of course, this objection can be raised for any 
perturbation calculation for strongly coupled linear sigma model). Even more importantly, 



as clearly stated in ref. ||22|, the Markov approximation (using which the dissipation term 



was derived) is valid in the weak coupling regime where the relaxation time of the soft modes 
is much larger than the collision time. In order that the soft modes have larger relaxation 
time, one can not take Vj as large as 40 fm~^ or so. 

As estimates of Vj from ref. p2|j8|] for the linear sigma model case are not reliable, it is 
reasonable to explore range of values of r/'. This is what we have done. For large values of Vj 
we get too much thermal fluctuations, while for smaller values thermal fluctuations are under 
control. In fact, one can argue that it virtually makes no sense to take Vj as large as 40 fm~^ in 
the Langevin equation. With such large 77', second order time derivative in the fleld equations 
becomes ineffective and the situation becomes that of extreme damping. This would not 
look like any relativistic system any more and would resemble more a damped condensed 
matter system. Also, in order that soft modes have at least somewhat larger relaxation time 
than the hard modes (which is the basis of whole analysis in ref. |^2|) it is important that r]' 
should not be much larger than 1 fm~^, so that soft modes do not thermalize immediately. 
For example when we express the equation in terms of dimensionless variables (so, the only 
parameter that enters fleld equations is A), we flnd that for our case, the value of r/' = 0.6 
fm~^ translates to a value about 0.5, whereas a value of 40 fm~^ will translate to about 32. 
This, we believe, also suggests that such large values of r( are not reasonable. 

With these points in mind, we go ahead with using Vj = 0.6 fm~^. We will also quote 
results for larger values of r]' (e.g. 1, 4 , and 10 fm~^) for comparison. We will see that 
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qualitative aspects of our scenario are not unaffected when we increase 77' to 4 fm~^. As 
thermal fluctuations are larger in this case, formation of very large DCC is less likely in this 
case. For example in one estimate we argue that for 77' ~ 4 fm~^ the probability of a large 
DCC is about 1/1000 (per event) compared to the case when 77' ~ 1 fm~^ 

Volume V over which the evolution of the soft modes is considered via Eqn.(3), will 
expand as plasma expands. We take V = Vo at r = tq. For the linear heating case, this 
amounts to taking a smaller value for V at r = rj„, which then expands to become equal to 
Vo at T = Tq. The volume expands via the following equation. 

It may appear reasonable to take Vo to be of order of a correlation volume. Then, a 
natural value to take will be Vo — 1 fm^. (This is for the instantaneous heating case. 
For the linear heating case it will lead to consideration of rather small volume at r = Tin, 
especially for the spherical expansion case.) Again, thermal fluctuations are larger in this 
case, so we give results for a larger value Vo = 5 fm'^. (The radius of such a region will be 
about 1 fm.) Results for even larger value of Vo further reduce thermally generated pionic 
components, otherwise results are essentially the same. It is important to remember that V 
does not really represent a correlation volume. It is a volume over which the average value 
of the chiral field is being evolved. Thus, in principle, we can take even larger values of Vo 
(see ref. 0, for example, where values of Vo as large as 10, 100, 1000 fm^ have been used). 
However, we are also interested in knowing the inhomogeneity in the chiral field generated 
due to thermal fluctuations during its overshooting to the opposite point of the vacuum 
manifold. In other words, we want to know the structure of the patch the chiral field forms 
near the south pole on in the entire region Q of the plasma after overshoot. For this, 
one needs to divide the entire region into several domains, and study the evolution of the 
average field in these domains. From this point of view, we should not choose V to be too 
large compared to the correlation length, as that will govern the scale over which thermal 
fluctuations will generate non-uniformities. For this reason, we let V expand only to a value 
of about 10 fm^, and freeze V to this value during further evolution. If we allow volume to 
keep expanding, or freeze it at a larger value, the effects of thermal fluctuations are further 
suppressed. 

We have used standard stochastic Runge-Kutta algorithm for solving the above Langevin 
equations which is second order in accuracy for the deterministic part of the equation and 



first order for the noise part |2^. There still remains the issue of the appropriate time step 
size for the noise term. For the case of white noise, as in our case, it seems clear that time 
step can not be arbitrarily small. We use time step size equal to 0.005 fm. We have checked 
how our results vary when this time step size is increased and decreased by a factor of 10. 
We find that the results with 77' = are unchanged, so the deterministic part of the equation 
is not much sensitive to the time step size, as one would expect (as long as the time step size 
is not too large to cause error in the evolution by the deterministic part). However, with 
r]' 7^ 0, the results depend more sensitively on the choice of time step size because of the 
noise term. Though, happily, we find that our main results do not change much, even when 
we vary time step size by such large amount. We find that the thermally generated pionic 
components, before the chiral field is de-stabilized by the appearance of the saddle point. 
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do not change much when we vary the time step size. Subsequent evolution of the chiral 
field varies significantly with such large variation of time step (we re-emphasize, for 77' = 
entire evolution had insignificant variation). However, this is easy to understand, as even 
a small change in the pionic components, near the saddle point, can lead to large variation 
during the subsequent rolling down of the chiral field on the vacuum manifold. In general, 
it seems that the choice of time step size can affect the results due to the assumption of 
white noise. In such situations, one will need to choose an appropriate value of the time 
step, presumably motivated by physics of the problem. For example, one may need to keep 
the time step size very small in order to get accurate evolution by the deterministic part of 
the equation. However, the noise term may be updated over a larger time step which could 
be motivated from the intrinsic scales in the problem. 

For the initial conditions, we choose the chiral field to be in the true vacuum at zero 
temperature, at r = 0, with the chiral angle 9 = 0, that is, $1,2,3 = and $4 = $i 
is taken to be zero initially. We will give plots for the case of spherical expansion of the 
plasma as thermal fiuctuations are too large for the case of linear expansion. For some cases 
we will also quote results for the linear expansion case. The evolution of a field is shown 
in Fig. 3 corresponding to temperature evolutions shown in Fig. 2. In two of the cases we 
see, in Fig. 3, that a field overshoots the central barrier and ends up becoming maximally 
disoriented with the chiral angle 6 being (close to) tt, (eventually rolling back to ^ = 0), 
e.g. the solid curve corresponding to the instantaneous heating case, with tq = 0.4 fm. For 
tq = 0.2 fm, (again, for instantaneous heating, with the same value of Tq = 350 MeV), the 
field does not overshoot as shown by the dashed curve in Fig. 3 and the field settles to ^ = 
without ever crossing the central barrier. This happens due to a faster cooling rate in this 
case, so that the force on the chiral field due to changing effective potential is operative only 
for a small time. Sometimes, the field settles down back at ^ = after crossing the central 
barrier more than once, as shown by the case corresponding to the linear heating case with 
Tin = 0.2 fm and tq = 0.4 fm (dotted curve). DCC formation can happen even in such 
cases, due to additional sources of dissipation in the field dynamics arising due to coupling 
to other fields, as well as due to particle production etc. We mention here that in our model 
Tq can, in principle, be below the critical temperature itself (see, ref. in this context). 
For example, for the instantaneous heating case with tq = 3.0 fm (which certainly is a rather 
large value) we find that the field overshoots to the opposite point on the vacuum manifold 
with To as small as 107 MeV. This we get for the case of longitudinal expansion, and in the 
absence of thermal fiuctuations, i.e. t]' = 0. With non-zero t]', thermal fiuctuations are too 
large for such a large value of tq. (For the spherical expansion case, we find the overshooting 
happens with Tq = 114 MeV for tq = 5 fm.) 

We emphasize here that our interest in a large DCC is from the point of view of its 
experimental signature in terms of measuring the ratio R. Thus it suffices that at every 
point in the physical space corresponding to a DCC, the chiral field rolls down along the 
same direction on S^, even if the field has some spatial/temporal variation. In fact, the 
experimental signature of DCC in terms of the distribution of the ratio R can not differen- 
tiate a DCC which has essentially uniform value of the chiral field, from a DCC where the 
chiral field has some spatial/temporal variation (still forming a condensate) but rolls down 
approximately along the same path on the vacuum manifold S^. One could distinguish the 
two cases if one could measure the coherence of pions coming out of the decay of the DCC. 
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For each set of parameters, we generate an ensemble of trajectories and find out the 
fraction of trajectories which overshoot the central bump at some temperature T < 115 
MeV (which is the point when the saddle point first appears). These are the trajectories 
which show how the chiral field gets disoriented from the true vacuum in the respective 
elementary domains in the region Q. Due to thermal fluctuations (with 77' ^ 0), pionic 
components become non-zero during the overshoot and the chiral field typically ends up at 
some point different from 9 = n. For the whole region to look like a single DCC, we require 
that the variation of the chiral field in the DCC arising due to thermal fluctuations remains 
small. We do not want to constrain pionic components which grow due to rolling down the 
valley after overshooting as that corresponds to the evolution of DCC. Thus we require that 
pionic components at the stage when the temperature is less than about 115 MeV, remain 
less than about 13 MeV (which is half of the shortest radius of the valley at that stage, i.e. 
the value of a at the saddle point with chiral angle equal to tt) , making sure that the average 
value of a for the trajectories overshoots the central bump. 

This cut off for the magnitude of pionic components is somewhat ad hoc. However, we 
would like to mention that the value of this cut off is only relevant in determining the allowed 
range of Tq. This range will be somewhat larger for a larger value of the cut off, but does 
not depend on it too sensitively. Fig. 4 shows parametric plot of the magnitude of 7? vs. 
(7, showing the evolution of the chiral field in the presence of thermal fluctuations. Here we 
have plotted different trajectories, difference arising due to the random nature of the noise 
term in Eqn.(3), corresponding to the situation with temperature evolution shown by the 
solid curve in Fig.2 (for which a evolution is shown by the solid curve in Fig. 3). As discussed 
earlier, these different trajectories should be thought of as representing the evolution of the 
(mean) chiral field in different portions, each of volume V, in the entire plasma region VL. 

This case represents the situation of a reasonably coherent DCC. Here, pionic components 
are dominated by thermal fluctuations upto a point where a ~ — 26Afey. This is also the 
stage around which the temperature takes the value of about 115 MeV, when the saddle 
point appears, making the chiral field unstable. At this stage one can say that DCC has 
formed. For later times, i.e. for a less than this value, the roll down of the chiral field starts 
dominating which is the situation of the evolution of the DCC. It is important to realize 
that much of the apparently erratic behavior of a given trajectory, in Fig.4, is not due to 
thermal fluctuations. It arises due to rolling down on the vacuum manifold, and due to 
oscillations about the valley during the course of roll down. Note that we plot a vs. |7f|. In 
this case we see that the average value of the pionic components due to thermal fluctuations 
at the formation stage of DCC, before the roll down of the chiral field starts dominating, 
remains less than about 13 MeV. This value is much smaller than the radius of the valley in 
the effective potential which has a value of about 26 MeV in = tt direction, and a value of 
about 51 MeV in = direction. More importantly, the area spanned by these thermally 
generated pionic components at this stage leads to a patch with area less than 3 % of the 
total surface area of the vacuum manifold at that stage, i.e. when the temperature is about 
115 MeV. 

Fig. 5 shows the scatter plot of the values of the pionic components at T = 115 MeV, when 
the saddle point at ^ = tt first appears. This scatter plot shows the projection of the values 
of the pionic components in the tti — 112 plane for 1000 trajectories. As we mentioned, the 
area of this patch in tti, 7r2, tts space is about 3 % of the surface area of the vacuum manifold 
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at that stage, which we approximate by taking the average of the areas of the 3-spheres 
with the smallest and the largest radii of the valley in the effective potential. This DCC 
then represents a large region in the physical space where the chiral field is (approximately) 
maximally disoriented, with the spatial variation of the chiral field in the entire region 
remaining relatively small. We mention here that evolution of individual regions with few 
fm size in Q via Eqn.(3), neglecting spatial variations, should be justified only up to r ~ few 
fm or so. For larger times, coupling between different domains should become important. 
For a proper, large time evolution, one needs to carry out a full simulation of the field 
evolution in Q. 

Thus we see that in our model, 6 distribution in the DCC typically represents a small 
patch S on the vacuum manifold S*^ at ^ = vr, i.e. near the south pole on S^. If this 
patch does not cover/enclose 6 = ir then the chiral field in the entire DCC should eventually 
roll down roughly along one direction (strip) on S^, leading to emission of (approximately) 
coherent pions, just as will be expected for any large DCC. Let us estimate the probability 
of forming such a patch near 6 = n whose evolution everywhere in Q leads to pions with 
roughly the same value of ratio R of neutral pions to all pions. Let us parameterize chiral 
field on by (a, tts, 112, vti) = {cos9, sinOcoscl), sinOsincpsinx, sin9sin(j)cosx)- By taking the 
number of pions to be proportional to the square of the corresponding field component, it 
is easy to show that R = cos^cf), see ref. For a rough estimate, divide range of (from 
to tt) in four equal parts. It is clear that if the chiral field in the patch S lies in any two 
parts, symmetric about = 7r/2, and its evolution remains confined to these same parts, 
then it should lead to roughly same average value of R for emitted pions. The probability 
for this type of patch S to form is roughly {\)^^~^\ where is the number of individual 
domains in Vt. Probability for other parameterizations of will be same. Since we freeze 
the volume V after initial expansion at a value of 10 fm^, the probability of forming a DCC 
with 100 fm^ volume will be about 0.002. For a 200 fm^ DCC this probability is about 10~^. 
For larger V probability for a large DCC is significantly higher. For example, if V is frozen 
after expansion to a value of 20 fm^, then the probability for a 200 fm^ DCC will be about 
0.002. It is possible to argue that similar evolution is expected to happen also when the 
patch covers 6 = it but is very asymmetric about it. Clearly, this can lead to a much larger 
DCC. A careful study of DCC evolution in our model, with appropriate boundary conditions 
at the surface of Q where chiral field is fixed at ^ = 0, requires numerical simulation of the 
chiral model in 3+1 dimensions, and we plan to present it in a future study. It is clear 
that generic situation will be when the patch is reasonably symmetric about 6 = it after the 
overshooting. This leads to some interesting possibilities and we will discuss it in detail in 
the next section. 

In the above estimate of the probability of DCC formation we did not account for the 
fact that due to thermal fluctuations the trajectories in some of the elementary domains in fl 
may not overshoot the central bump (after the time the temperature drops below 115 MeV). 
We have checked that with with t]' = 0.6 fm~^, essentially all (99%) trajectories overshoot 
the central barrier after temperature drops below 115 MeV, so there is no extra suppression 
in the probability of DCC formation. With rj' = 1 fm~^ also, almost all (97 %) trajectories 
overshoot. When we repeated the simulation for a larger value of t]' (= 4 fm~^), we find 
that thermal fluctuations are larger and a smaller fraction (~ 45%) of trajectories overshoot 
the central bump. (In all these cases we keep Tin = 0.4 fm and T^ax = 350 MeV, with 
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instantaneous heating). This will suppress the probability of formation of large DCC For 
example, in the above estimate of the probability, we should multiply by a factor of about 
1/2 for each elementary domain (in the region Q) to account for the probability that the 
trajectory in that domain may not overshoot the central barrier. This means that for a DCC 
of volume 100 fm^, with each elementary domain having 10 fm^ volume, the probability will 
be suppressed by an additional factor of 10~^. Also, due to larger thermal fluctuations in this 
case, at the time trajectories overshoot the central barrier, the average pionic components 
form a patch of area which is about 27 % of the area of the vacuum manifold. This represents 
a DCC which is not very coherent, though its subsequent evolution may still be similar to 
a DCC. 

When we increase t]' to a value as large as 10 fm~^ (which is about 1/4 of the estimate 
based on ref. fS^J^), we find that the patch formed by pionic components (at the time of 
overshooting the central bump) is too large, covering almost half of the vacuum manifold 
and destroying the coherence of DCC. Also, now a very small fraction, about 8 %, of the 
trajectories overshoot the central barrier after temperature drops to 115 MeV which will 
lead to strong suppression in the probability of formation of large DCC. 

We have also checked how these results change if we take different set of parameters. 
For example, for the set of parameters used in ref. (with A = 20), we find that Tin = 0.9 
fm and Tmax = 280 MeV is preferred. With these, about 40 % trajectories overshoot the 
central barrier for rj' = 0.4 fm~^, with pionic components at overshoot forming a patch of 
area which is about 15 % of the area of the vacuum manifold. For larger rj', the fraction of 
trajectories overshooting becomes smaller (about 20 % for rj' = 1 fm~^). 

Having discussed the effects of larger values of r]', let us now discuss constraints on the 
value of To in our model. For small Tq, the chiral field does not pick up enough energy to 
overcome the central barrier. We find the lowest allowed value T^m of Tq to correspond to 
the case of linear heating with tq = 0.6 fm, and rj„ = 0.3 fm. For this case T^m is found 
to be about 225 MeV. With tq = 0.5 fm we find T^m to be about 240 MeV for the linear 
heating case, (for the instantaneous heating case with tq = 0.5 fm we find T^m = 260 MeV). 
Tjnin increases with decreasing value of tq. For values of tq larger than 0.6 fm, thermal 
fluctuations become dominant and there is no allowed range for Tq. 

Upper bound T^ax on Tq arises from the requirement that thermal fluctuations do not 
dominate, so the entire region still looks like an approximate DCC. For this we impose re- 
quirement that thermally generated pionic components remain relatively small. As discussed 
earlier, we do not want to constrain pionic components which grow due to rolling down the 
valley after overshooting. Thus we require that pionic components at the stage when the 
temperature is 115 MeV, remain less than about half of the shortest radius of the valley at 
that stage which is about 26 MeV. Largest value of Tmax is found to be about 390 MeV, for 
the case of instantaneous heating with tq = 0.4 fm. With tq = 0.5 fm, Tmax is obtained to be 
about 310 MeV (again, for the instantaneous heating case). Tmax decreases with increasing 
value of Tq. Values of tq smaller than 0.4 fm appear not very realistic, hence we do not 
consider such small values. Thus, flnally we get the allowed range for Tq to be roughly 225 - 
390 MeV. The plots in Fig.4 and Fig.5 have been given for the case of instantaneous heating 
case with tq = 0.4 fm. The allowed range of Tq for that particular case is about 290 - 390 
MeV. 

These values of Tmax are significantly lower than the value of the QGP temperature 
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expected in central events at heavy-ion collisions at LHC and RHIC involving heavy nuclei. 
This suggests that in order to look for spectacular signals of DCC in these experiments one 
should analyze events with a range of values of impact parameter, i.e. lower multiplicity 
events, as the QGP temperature for large impact parameter collision will be expected to 
be smaller. Smaller nuclei should be more appropriate from this point of view. In fact, it 
seems that AGS and SPS may be better suited for large DCC, especially if one analyzes 
events with a range of centrality. In this way some events may have Tq falling in the right 
range for DCC formation. However, the requirement of early spherical expansion will prefer 
p — p {p — p) collisions for our scenario to be effective. In this sense, Tevatron may indeed 
be a good place to look for DCC 0]. This seems entirely consistent with the observation 
of Centauro events where heavy nuclei are not expected to be involved. We may mention 
here that even for the conventional models of DCC formation, where DCC domain sizes 
are expected to be small, it seems better to probe peripheral events. This is because for 
peripheral events, resulting QGP region will be expected to have small volume. This will 
allow very few DCC domains to form in a single event, hopefully making it easier to analyze 
signals of DCC pions. Importance of peripheral events has also been emphasized in ref. |]2^ 
based on somewhat different considerations. 



IV. TOPOLOGY OF FIELD CONFIGURATION AND PION DISTRIBUTION 

Let us now discuss the case when the patch formed near the south pole on S^, after 
overshooting of the chiral field, is reasonably symmetric. As we mentioned earlier, this 
type of situation will be expected to be most generic. It is important to note that even 
when thermal fluctuations are much more dominant than considered above, this type of 
field configuration will still form. This is essentially because larger thermal fluctuations will 
only increase the size of the patch at the south pole (for example in the case of 77' = 10 fm~^, 
as discussed above). This is very important, especially in view of the uncertainties in the 
estimates of rj' (as discussed in the previous section). As we showed earlier, for large values 
of rj', formation of a reasonably coherent, large DCC becomes less likely. However, even in 
that case, our considerations of this section will apply. 

When the patch near the south pole is reasonably symmetric (after overshooting) then 
we expect in this case that the field will roll back towards the north pole, which is the true 
vacuum, in different directions. This will lead to covering of entire 5*^ once, and hence, will 
lead to the formation of a single large Skyrmion like configuration (due to fixed value ^ = 
outside fl, see ref. ||26|). As such it does not look anywhere hke the standard picture of 
evolution of a large DCC. However, we now argue for the possibility that the evolution of 
such a configuration leads to spatial distribution of pions which can be grouped together in 
terms of definite values of R, thereby directly probing disoriented condensates of the chiral 
field. With the parameterization of chiral field used above, R depends only on angle (f) on 
the vacuum manifold {R = cos^(0)). 

Now, given the situation that we have a (very) large Skyrmion like configuration in the 
region Q (possibly of size 15 fm or so), the configuration should evolve to a spherically 
symmetric configuration of standard Skyrmion |^ during the process of its collapse, for 



energetic reasons. This reasoning is supported by numerical simulations of texture formation 
and collapse in the early Universe Textures are just like Skyrmions, apart from the 
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effect of pion mass term. One can argue that the pion mass term should not affect these 
considerations of texture (Skyrmion) collapse. In ref. [27|, it was found that a texture of 
generic initial configuration becomes spherically symmetric, and approaches the analytical 
solution, during collapse. 

Let us assume the region Q or, more precisely, a smaller region containing the collapsing 
Skyrmion like configuration, to be spherical (which should be fairly justified for late stages 
of evolution) with a and (3 being the polar and azimuthal angles on the surface of Q (which 
is a two sphere 5*^). Let us further assume that the Skyrmion configuration corresponds 
to angles 6, 0, x on the chiral vacuum manifold 5*^ varying spatially so that (p = a, x = P 
and 6 varying radially outward from 6 = n at the center of fi, to 6* = outside Q, which 
is the true vacuum. With such a Skyrmion configuration, we see that the chiral field in a 
conical shell in fl, corresponding to constant value of angle a, has fixed value of and hence 
corresponds to the same value of the ratio R. 

If we assume that DCC pions will come out predominantly radially from Q due to rapid 
expansion of the region, which does not seem unreasonable as these pions should also have 
very small momentum in the rest frame of the plasma, then we simply need to collect pions at 
the detector with fixed polar angle a. Since one does not know, a priori, about the orientation 
of the Skyrmion configuration in Q, one should analyze pions in conical shells with fixed a 
and with varying choices of the north pole of the two sphere in a Air detector around the 
collision point. It is easy to see that this will also take care of other parameterizations of the 
chiral field which are different from the one used in the above argument. For one particular 
choice of the north pole one will start getting pions with R varying as cos'^a, for each choice 
of circle on S^, say, in a small interval Aa at fixed a. This analysis can be suitably adopted 
for fixed target experiments. We would like to mention that this technique of analyzing 
pion distribution seems interesting by itself. If there is any axially symmetric phenomenon 
happening in heavy-ion collisions which leaves imprints on the distribution of low pt pions, 
then this analysis may be able to detect it. 



V. CONCLUSIONS 

We summarize our main results. We have considered the dynamics of chiral field dur- 
ing very early stages of the evolution of the plasma. We model the early, pre-equilibrium, 
evolution of the system in terms of finite temperature effective potential, with the temper- 
ature first rising to a maximum value Tq (the heating stage), and then decreasing due to 
continued plasma expansion. Due to change in the effective potential during the heating 
stage, the chiral field is pushed over to the opposite point on the vacuum manifold, thereby 
maximally disorienting the chiral field in the entire plasma region. This should convert the 
entire plasma region into a single DCC domain. Thermal effects tend to create fluctuations 
in this DCC region. The requirement that such fluctuations do not become dominant, so 
that the whole region still approximates a single DCC domain, puts constraints on the rate 
of plasma expansion, as well as on the maximum allowed value of Tq. We find that spheri- 
cal expansion is preferred over the longitudinal expansion. During very early stages of the 
expansion, for times much less than the transverse dimension, it does not seem reasonable 
to invoke spherical expansion. Thus for heavy nuclei, formation of such large DCC domain 
looks less likely. From this point of view, hadron-hadron collisions should be a good place 
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to look for the formation of large DCC since there spherical expansion will be expected to 
dominate from rather early stages. 

Even for the spherical expansion case, if Tq is too large then the system spends too 
much time at high temperatures, again leading to large thermal fluctuations. We find that 
thermal fluctuations remain relatively small only when Tq remains less than about 390 MeV. 
There is also a lower bound on Tq arising due to the fact that for small Tq the chiral field 
does not overshoot at all. We estimate this lower bound to be about 225 MeV. There are 
uncertainties in the determination of these bounds, for example due to the value of cutoff of 
pionic components at overshoot etc. Conservatively, one can say that there is a window for 
Tq from about 200 MeV to about 400 MeV, which allows for the formation of such a large 
DCC. We emphasize again that these values are rough estimates. Evolution of the field in 
one region of the plasma will be affected by the evolution of the field in the neighboring 
regions. This requires a full 3-dimensional simulation, with proper boundary conditions. 

We have also proposed an interesting possibility for the field evolution in the entire 
plasma region. Here the field, after overshooting, evolves to form a topologically non-trivial 
structure which looks like a large Skyrmion. We then argue that the collapse of this structure 
should lead to a very definite type of spatial distribution of pions. By proper grouping of 
pions, this should lead to definite values of the ratio R arising in a predictable manner 
in a single event. We suggest that this technique of analyzing particle distributions could 
be fruitful, not just in the context of our model, but in uncovering any axially symmetric 
phenomenon in these experiments which leaves imprints on particle multiplicities. 

We mention that the ideas discussed here have strong implications for condensed matter 



systems. One example is the formation of defects via a new mechanism, discussed in ref. p8 
where field oscillations and subsequent flipping (which is just like the chiral field overshooting 
to the opposite point on the vacuum manifold) leads to defect formation. Thus, if a small 
portion of superfluid helium, or superconductor, is heated rapidly to a temperature less 
than the critical temperature, then it may lead to topological defect formation, without ever 
going through a phase transition. We will present these results in a separate work. 
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FIGURE CAPTIONS 



1) Thick solid curve represents the effective potential at T = with the value of the 
chiral field denoted by point p. Temperature is instantaneously raised to T'i(< Tc). Dashed 
curve shows the effective potential at this temperature. 

2) Evolution of the temperature T (in MeV) as a function of proper time r (in fm). 
Dotted curve shows hnear heating. Sohd and dashed curves show instantaneous heating 
with different values of tq. 

3) The evolution of a (in MeV). r is in fm. 

4) Different trajectories of the chiral field for the instantaneous heating case (solid curve 
in Fig. 2). Figure shows parametric plots (upto r = 10 fm) of a and |7f| (both in MeV). 

5) Scatter plot of the projection of pionic components in the tti —ti2 plane corresponding 
to Fig.4, at the stage when T = 115 MeV. This patch covers only about 3 % of the area of 
the vacuum manifold at this stage. 
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FIG. 2. 
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FIG. 3. 
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FIG. 4. 
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FIG. 5. 
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